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Preface 


This  is  the  suminary  report  for  Al-'OSR  grant  76-3109  "N'onparamctric  Methods 
with  Applications  to  Reliability."  The  grant  period  was  October  1,  1976 
through  September  30,  1978.  [Beginning  October  1,  1 7 , the  principal 
investigator  Myles  Hollander  continues  research  as  co-principal  investigator 
(with  Frank  Proschan)  under  AFOSR  grant  76-3678  "Statistical  Aspects  of 
Reliability,  Maintainability,  and  Availability."] 

In  tli is  summary  we: 

a.  List  specifically  the  accomplishments  in  terms  of  technical  reports 
that  have  been  written,  and  articles  that  have  been  published  or  are  scheduled 
to  appear. 

b.  Give,  for  each  paper  and/or  technical  report,  a technical  summary 
(abstract),  anu  a non-tochnical  summary  with  the  latter  describing  potential 
applications  to  the  Air  Force. 

c.  List  the  specific  activities  and  honors  of  note  of  the  principal 
investigator  during  the  contract  period. 


1.  List  of  AccompI  i shmeiits 
Technical  Reports  Produced  Under  Grant  AIOSR  76-3109 

1.  Nonparamet ric  empirical  Hayes  estimation  of  the  probability  that  X V, 
by  N'.  Ho  1 lander  and  It.  M.  Korwar,  #1,  October,  1976. 

2.  'lest  inn  for  agreement  between  two  groups  of  judges,  by  M.  Hollander  and 
,1 . Sethuvainan,  #2,  .January,  1977. 

3.  Nonpara met ric  Bayes  estimation  with  incomplete  Dirichlet  prior  informat ion, 
by  (I.  Campbell  and  M.  Hollander,  f 3,  June,  1977. 

I.  Prediction  intervals  with  the  Dirichlet  prior,  by  C.  Campbell  and  M. 
Hollander,  M,  June,  1977. 

3.  Testing  whether  more  failures  occur  later,  by  M.  Hollander,  #f>,  October, 

19/7. 

(>.  Testing  to  determine  the  underlying  distribution  using  randomly  censored 
data,  by  M.  Hollander  and  1;.  Proschan,  #6,  May,  1978. 

7.  On  the  normal  convergence  of  a family  of  simple  epidemics,  by  \.  A.  l.angberg, 
#7,  June,  1978. 

8.  The  discrete  asymptotic  behaviour  of  a simple  batch  epidemic  process, 
by  L.  Billard,  H.  l.acayo,  and  N.  l.angberg,  I1 8,  June,  1978. 

9.  On  the  negative  binomial  convergence  in  a class  ol  m-dimensional  epidemics, 
by  11.  l.acayo  and  N.  L.angberg,  1/9,  July,  1978. 

10.  The  symmetric  m-dimensional  simple  epidemic  process,  by  L.  Billard, 

II.  l.acayo,  and  N.  l.angberg,  MO,  June,  1978. 

11.  Characterization  of  nonparametric  classes  of  life  distributions,  b\ 

N.  A.  l.angberg,  II.  I.cfin,  and  F.  Proschan,  #11,  July,  1978. 
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1-’.  On  .1  cha ract er i at  ion  ot  mill t i var  i at c distributions  i,  11  h applications 
in  reliability  hiuI  epidemiology,  by  N.  A.  l.angberg,  ’1',  Align  - 1 , 197S. 

Articles  and  Hook-'.  (Published  or  In  I’ross) 
l.  Nonparamot  r i >.  empirical  Bayes  estimation  ofthe  prnbubi  1 i ty  that  X Y, 
by  M.  Ilollamlor  and  K.  M.  korwar.  Co  nr  mui  i cat  ions  in  Statistics  Tbcor. 

Moth.,  AS 1 14)  , 1369-1383,  1976. 

Statistics:  A B iomedical  Introduct ion  by  B.  Iv.  Brown,  Jr.  and  M.  Hollander 
A volume  in  the  Wiley  Scries  in  Probability  and  Statistics,  John  Wiley,  New 
York,  Hi 77. 

3.  Pune t ions  decreasing  in  transposition  and  then  applications  in  ranking 
problems,  by  M.  Hollander,  I.  Proschan,  and  d.  Sethuraman . Annals  of  Statistic 
722-733,  1977. 

•I . Rank  order  estimation  with  the  Hirichlet  prior,  by  G.  Campbell  and  M. 
Hollander.  Annals  of  Statistics  (> , 112-153,  11)78. 

5*.  Nonpar. nnet ric  estimation  of  distribution  functions,  by  M.  Hollander  and 
k.  M.  Korwar . The  Theory  and  Applications  of  Reliability,  Vol . ed.  by 
C.  P.  fsokos  and  l.  N.  Shimi,  Vc.ademic  Press,  New  York,  85-107,  1977. 

().  resting  whether  more  failures  occur  later,  by  M.  Hollander.  Proceedings 
of  the  1978  Annual  Reliability  and  Maintainability  Sympos ium,  103  I (>t* , 1978. 

7.  resting  for  agreement  between  two  groups  of  judges,  by  M.  Hollander  and 
.1 . Sethuraman.  To  appear  in  Bioroetriki  65,  1978. 

8.  Nonpa rametr i c Hayes  estimation  with  incomplete  Hirichlet  prior  information, 
by  G.  Campbell  and  M.  Hollander.  To  appear  in  Optimising  Methods  in  Statistics 
cd.  by  .1.  Rustagi,  Academic  Press,  New  York,  1979. 

* This  paper  was  not  written  or  revised  during  the  contract  period  and  thus  it 
does  not  acknowledge  Grant  /<>  3109.  It  is  included  here  for  completeness 
since  it  is  a publication  of  the  principal  investigator  which  appeared  duriny 
the  contract  period. 


2.  Pi scuss ion  of  pa pc rs . 

. Nonparan  trie  empirical  Bayes  estimation  of  the  probabi lit)  tj  it  X • V. 

2 . 1 . a . Abst ract . 

let  X and  V be  two  real  valued  independent  random  variables  with  distribution 
tune!  ions  I-  and  (i,  respect  ivel>  . Ke  consider  the  problem  of  est  innt  ins’,  tin. 
probability  that  X S V,  denoted  by  A,  A - / IMG.  Motivated  by  I'erguson's  (1973) 
nonpa rame t r i c Bayes  estimator,  we  propose  an  empirical  Bayes  estimator  of  A which 
requires  less  prior  information  about  ot  ^ (_  * ) and  x (•},  the  parameters  of  the 
(assumed)  underlying  Pirichlet  processes.  Onl)  UjfK)  and  u,(!')  need  be  specified. 
Consider,  then,  the  following  set  up  appropriate  for  an  empirical  Bayes  estimation 


problem.  Let  {xP1},  lxf2)},  i = 1,2, 

‘ t i 


be  two  independent  sequences  of  iiule- 
c i ,,  in,,..,  vC.n  . / .(j) 


pendent  vectors  of  observations  from  )■  and  (i  respectively.  Here  X;-  - ( X:  - ‘ , . . . , X^  ) , 

"'j 

i 1,2  and  i = 1,2,...  . Assume  independent  Pirichlet  priors  on  (w,  K ) with 
parameters  and  respectively  for  I and  C.  Here  R is  the  real  line  and  (5 
is  the  o field  of  Bore]  subsets  of  R.  l.et  the  action  space  be  the  closed  interval 

A ~>  „ 

[0,  1 | , and  the  loss  function  be  I . ( A , A)  (A  - A)^,  where  A is  an  estimator  of  A. 

We  assume  (E)  and  u -,  ( R)  arc  known.  We  then  propose  the  estimator  ' be  low  as 
an  estimator  of  A on  the  (n+l)tli  occasion.  The  estimator  is  given  by 


A = A(X[1) xt1];  X[2> X^j) 

n I n+1  1 ’ n+l 


n n , n n, 


P,  n P2  n ~2~  l l l l l ,2 

’ 1 ’2  n n.n  , i-1  j = l k I f--l  ( X>  T 

I / * V I 1 


(XlU) 
)/  kf 


+ 1’l.,.1(Up2.nJ)  .uyf,  J,  ^,%U)((_ 


1 I 1 


J I lk 


Ik 


11 1 n2 

Cl  p , Hi  - |>,  ) 1 ) } l 

' -»'H  > ; -I  1 ! 


1 ' I j I ( -W,X 


lit  I .1 


ix,ll  ) 

n * 1 , » 


where  p.  ^ u.  lR)/lo  . (k‘)  ♦ in),  i 1,2,  uiul 

1 if  x t A 

A (A)  < 

0 if  \ { A. 

2.1. b.  Potential  use  to  Air  Force. 

Fstimntion  of  P(.\  < X ) , the  probability  that  X is  less  tliat  V,  is  part  icu 
larly  important  in  engineering  situations  where  X is  stress  ;nul  Y is  strength. 
Applications  ahouiul  in  materials  development  and  evaluation  situations.  for 
example,  are  particular  noset ips  and  heat  shields  of  reent r>  vehicles  sufficiently 
strong  to  withstand  erosion  due  to  impact  of  rain,  ice.  debris,  etc.?  Here  the 
stress  is  created  partial!'  by  the  thermal  environment.  fhe  estimator  given 
by  Cl)  enables  effective  utilisation  of  available  data  (present  and  past  samples') 
without  the  imposition  of  overly  restrictive  distributional  assumpt  ions. 

2.1. C.  Status  of  paper . 

This  has  appeared  as  an  inv  ited  paper  in  a special  issue  liieceaiber,  IP  '«)  , 
of  lloi. mum  i cat  ions  in  Statistics,  devot  ed  to  recent  advances  in  nonpnramet  ric 
met  liv'd s . 

2.2.  Testing  for  agreement  between  two  groups  of  judges. 

2 . 2 .  a . Ah st  rac t . 

Ihe  "problem  of  in  rankings",  so  named  by  Kendall  and  studied  extensively 
by  Kendall  and  Rubington  Smith  ( IPoP) , Kendall  (.1PY0),  ami  others,  considers  the 
relationship  between  the  rankings  that  a group  of  ; nidges  assigns  to  a set  of  1 


objects.  In  this  paper  we  suppose  there  are  two  groups  of  in  ranking  l lie 

objects,  nitil  we  addro  . the  following  question  (posed  to  us  b>  Kendal  1J:  liiven 
that  there  is  agreement  within  each  group  oF  judges,  how  can  one  test  for  evidence 
of  agreement  between  the  two  groups?  A test  of  agreement  proposed  by  Schuc  in> 
and  f raw  ley  (1973),  and  further  advanced  hy  1 i and  Schucany  (1973),  is  shown  to 
be  misleading  and  does  not  proiidc  a satisfactory  answer  to  k-ndall's  question, 
he  then  adapt  a procedure,  proposed  hy  h'ald  and  Kolfowit.-.  (19-il)  in  a slightly 
different  context,  to  furnish  a new  test  for  agreement  between  two  groups  of 
j udges . 

2.2. t>.  Potential  use  to  Air  Force. 

A team  of  m experts  at  Aerospace  Laboratory  1 ranks  k 4 escape  systems. 

Hank  1 is  assigned  to  the  most  preferred  escape  system,  rank  to  the  second  pre- 
ferred, and  so  forth,  so  that  rank  4 is  assigned  to  the  least  preferred.  A 
different  team  of  experts  ranks  the  same  4 escape  systems  at  Aerospace  Lab  2. 

It  is  found  that  there  is  general  agreement  • each  team,  and  it  is  then 

desired  to  test  for  agreement  between  teams.  This  can  be  done  by  the  distribution 
free  test  we  have  developed. 

2.2. c.  Status  of  paper. 

This  paper  is  scheduled  to  appear  in  B i omet r i ka , 197S.  Sir  Maurice  Kendall, 
in  private  correspondence,  lias  termed  the  paper  a major  breakthrough. 

2.3.  Nonparametr ic  Bayes  estimation  with  incomplete  Dirichlct  prior 
in  format  ion . 

2. 3.  a.  Abstract . 

In  this  paper  we  treat  the  topic  of  incomplete  information  regarding  the 


J 


parameter  at  of  a birieb’et  process  prior.  Ferguson  (1973)  introduced  the  Dirichlct 
process  for  the  incorporation  of  prior  information  into  the  anal > sis  of  nonpar. imet ri 
problems.  The  process  can  be  viewed  as  a prior  on  the  set  of  all  distributions 
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i ■ 


i 

\\ 
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on  a meusinab  lo  space  (X,  Aj  . The  process  is  paramet  i i zed  liy  a,  a non  nyativ*--, 
non-null  finite  measure  on  (X,  A).  (In  this  paper  we  restrict  to  situations 
where  X = R,  the  real  line,  and  A B,  t lie  Borel  o-field.)  Typically,  to  use 
estimators  which  arc  Bayes  with  respect  to  a Dirichlet  process  with  parameter  a, 
the  statistician  must  provide  a complete  specification  of  the  measure  a.  This 
paper  develops  some  estimators  that  rely  only  on  parti;  1 information  concerning 

One  approach  to  incomplete  information  concerning  a is  that  initiated  by 
DoKsum  (11)72).  Do k sum  assumes  that  a(L-,  t . + ^ ) . i - l,...,k-l  are  Known  with 
a(R-(tj,  i | ) = 0.  That  is,  t lie  values  that  assigns  to  the  h-J  intervals 

(tj.t,] (tj  I’tpl  :irc  known,  and  u(R)  = a(tj,t^].  In  this  paper,  Uoksuin's 

technique  for  obtaining  a mixed  rule  is  considered  and  shown  also  to  yield  a 
C-minimax  rule  for  a suitable  choice  of  (!.  The  •approach  is  then  applied  to 
obtain  a new  mixed  rule  for  estimating  Pr(X  < Y)  . 

2.A.b.  Potential  use  t.o  Air  Force. 

The  general  approach  allows  one  to  incorporate  prior  information  into  the 
decision  problem.  Here  the  prior  ii  irmation  is  incomplete.  In  the  case  of 
complete  prior  information,  the  prior  guess  rt(*)/a(Rj  of  the  underlying  distri- 
bution is  specified  completely  whereas  in  this  approach  the  scientist  need  only 
"guess"  or  specify  a on  intervals.  The  noset ips  example,  discussed  in  2.1.1), 
is  applicable  here  since  we  have  derived  a mixed  estimator  for  Pr(\  s Y) . 
However,  the  mixed  estimator  only  requires  partial  prior  information. 

2.3. c.  Status  of  paper. 

This  paper  will  appear  in  Optimizing  Methods  in  Statistics  - 11.  ed.  by 
J.  Rustagi,  Academic  Press,  New  York,  1979. 

2.4.  Predict  ion  intervals  with  the  Dirichlet  prior. 

2.4.  a.  Abstract ■ 

I.et  Xj,...,X  and  Yj,...,Y  be  consecutive  samples  from  a Dirichlet  process 
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- 8 - 

on  (.R,  8)  (tin-  real  line  R with  the  Borel  c field  i>)  with  parameter  1 ypj  cal  1/ , 
prediction  intervals  employ  the  previous  observations  X ,X  iri  order  to  pro- 

1 • 1 ii 

diet  a specified  function  of  the  future  sample  Here  one  and  two  ided 

prediction  intervals  for  at  least  1.  of  N future  observations  are  developed  for 
the  situation  in  which,  in  addition  to  the  previous  sample,  there  is  prior  infor- 
mation available.  The  information  i s specified  via  the  parameter  a of  tin  Dirichlet 
process . 

2.4. h.  Potent  ial  use  to  Air  l-orce. 

In  many  practical  problems  in  quality  control  and  statistics,  one  wishes  to 
use  the  results  of  a past  sample  to  construct  an  interval  which  will  contain  the 
results  of  a future  sample  from  the  same  population  with  a specified  probability. 
Such  an  interval  is  known  as  a prediction  interval.  Prediction  intervals  are  of 
special  interest  to  engineers  who  are  concerned  with  setting  limits  on  the 
performance  of  a small  number  of  units  of  a product.  Such  limits  would  he 
required,  for  example,  in  setting  specification  to  contain  with  a high  proba- 
bility a critical  performance  characteristic  for  all  units  in  an  order  of  three 
heavy  transformers  when  the  only  available  information  is  the  data  on  five 
previous  transformers  of  the  same  type.  By  using  the  limits  of  a prediction 
interval  as  specification  limits,  one  is  assured  with  a desired  probability  that 
all  three  transformers  will  meet  specifications. 

Typically,  prediction  intervals  are  parametric,  Bayesian  parametric,  or 
(non-Bayesian)  nonparametric.  The  prediction  intervals  we  derive  are  Bayesian 
nonparametric,  and  thus  simultaneously  allow  the  user  to  incorporate  prior 
infromation  into  the  procedure,  without  requiring  restrictive  assumptions  about 
the  form  of  the  underlying  population. 

2 . 4 .  c . St atus  of  paper. 

The  paper  has  been  submitted  to  the  Journal  of  the  American  Statistical 
Assoc iat ion . It  was  well-received,  subject  to  minor  suggested  revisions  which 


arc  currently  being  made. 
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2.  S.  Functions  decrea  ting  in  transposition  and  their  appl ical  ions  in 
ranking  prohl cros. 

2. 5.  a.  Abstract . 

l.et  A = f A A ),  A,  s ...  < X , and  x = (x x ).  A function  ‘.(A,  x) 

is  said  to  be  decreasing  in  transposition  (!)T)  if  (ij  g is  unchanged  when  the  same 
permutation  is  applied  to  A and  to  x,  and  (iij  gCA,  x)  ? g(A,  x’)  whenever  x’ 
and  x differ  in  two  coordinates  only,  say  i and  j,  fx.  - x.)(i  - j)  > 0,  and 
x!  = x.,  x!  = x..  The  111'  class  of  functions  includes  as  special  cases  other  well 
known  classes  of  functions  such  as  Schur  functions,  totally  positive  functions  of 
order  two,  and  positive  set  functions,  all  of  which  are  useful  in  many  areas  in- 
cluding stochastic  comparisons.  Many  well  known  multivariate  densities  have  the 
DT  property.  This  paper  develops  many  of  the  basic  properties  of  DT  functions, 
derives  their  preservation  properties  under  mixtures,  compositions,  integral 
transformations,  etc.  A number  of  applications  are  then  made  to  problems  involving 
rank  statistics. 

2.5. b.  Potential  use  to  Air  Force. 

The  standard  assumption  in  comparing  failure  times  of  two  systems,  X,,...,X 

1 m 

of  system  1,  and  Y^,...,Yn  of  system  2,  is  that  the  X's  are  a random  sample  from 
one  population,  that  the  Y's  are  a random  sample  from  a second  population,  and 
that  all  X's  and  Y's  arc  mutually  independent.  Then  to  test,  for  example,  the 
hypothesis  of  equivalent  populations  against  the  alternative  that  one  system  tends 
to  have  shorter  failure  times  than  the  other,  two  sample  rank  statistics  can  he 
used.  However  in  many  Air  Force  applications,  the  assumption  of  independence 
can  be  violated,  say  by  having  the  systems  operating  in  a common  environment. 

The  DT  paper  allows  for  application  of  these  rank  tests  even  in  the  case  of  such 
dependency  situations  by  providing  inequalities  on  power  functions. 


I 0 

2.5.c.  St . * 1 1 1 * . of’  paper . 

lhe  paper  appeared  m I lie  July,  1077  issue  of  th  \tin.t  1 ol  Slat  i t ii  s, 

2.(>.  Rank  order  estimation  with  the  Oirichlct  prior. 

2 . f< . a . Abst  ract  . 

Suppose  that  a sample  of  si  e it  from  a J i str  i hut  ion  function  I is  obtained. 
However,  only  rl<  n)  values  from  the  sample  are  observe*!,  say  \ ,...,X  . Kithout 
loss  of  generality,  we  van  consiiler  X ,...,X  to  be  the  first  i values  in  the 
(unonlereJ)  sample.  I'he  problem  is  to  estimate  the  rank  order  (i  of  among 

X . The  situations  of  interest  include  f non-random,  either  known  or 

1 n 

unknown , and  f random.  the  random  case  assumes  that  I is  a random  distrihut ion 
function  chosen  according  to  ferguson's  (197 3)  l'irichlet  process  prior.  Since 
this  random  distribution  function  is  discrete  with  probability  one,  average  ranks 
are  used  to  resolve  ties.  A Hayes  estimator  (sipuired- error  loss!  of  (1  is  do  - 
ve loped  for  the  random  model.  for  the  non  random  distribution  function  model, 
optimal  non  Hayesian  estimators  are  developed  in  both  the  case  where  I is  known 
and  the  case  where  1 is  unknown.  Those  estimators  arc  compared  with  the  Hirichlct 
estimator  on  the  basis  of  average  mean  square  errors  under  both  the  random  and 
non- random  models. 

2.(>.l>.  Potential  use  to  Air  force. 

Consider  for  example,  the  following  situations: 

(.  i ) An  astronaut  ( IV  Iv , say)  undergoes,  as  one  of  a pilot  group  of  IS 
astronaut  trainees,  extensive  preparation  for  a space  mission.  I ach  asl  ivnaut 
earns  a score  X,  a measure  of  overall  performance.  Will's  score  is  \ . based  on 
the  observed  values  X|,...,Xj,,  we  wish  to  estimate  WW's  rank  in  the  total  pool 
of  50  trainees.  (Unlj  the  best  ten  astronauts,  as  measured  by  X,  will  be  chosen 


for  t h • mi ss i on . ) Mere  r 


15  and  n 50. 


I i i ) The  M.mt  ilia  Itiver  lias  II  ood  Oil  four  1 hues  in  this  decade  wit  It  t 1 1 
severity  i'f  each  llmul  measured  by  X,  l lu‘  height  of  l lie  river.  On  the  basis 

of  the  observations  X ,...,X|,  how  can  wo  estimate  the  severity  of  the  fiist 

flood,  in  the  group  of  these  four  and  the  next  five  that  occur?  Cquivalently, 

how  can  we  estimate  the  rank  order  of  Xj  among  \j \t)?  (Note  that  wc  could, 

for  example,  interchange  the  roles  of  X^  and  X,  and  pose  the  quest  ion  in  terms 
of  estimating  the  severity  of  the  fourth  flood.)  Here  r 4 and  n --  9. 

Of  course  example  (ii)  is  easily  generalised  to  cover  other  undesirable 
(or  desirable)  events,  and  example  (i)  is  applicable  in  other  situations  where 
a subgroup  is  selected  to  he  on  a team  or  perform  a mission.  The  rank  order 

estimators  derived  can  be  used  in  these  prediction  problems. 

2.6. C.  Status  of  jiajter. 

The  paper  appeared  in  the  January,  1978  issue  of  the  Annals  of  Statistics. 

2.7.  Nonparamet r ic  estimation  of  distribution  functions. 

2. 7.  a.  Abstract. 

Consider  the  nonpa ramet r ic  estimation  of  n + I distribution  functions 

F,,F, F .,  on  the  basis  of  samples  X.  (X.......X  ) from  1'.,  i I,. ...ml 

12  ’ n *■  1 1 t 1 1 i m . i 

In  an  empirical  Haves  context  where  X X were  the  "past"  samples,  korwar  and 

Hollander  ( 1D7(>)  used  Ferguson's  (197S)  Piriehlet  proeess  to  propose  an  estimator 
of  !■  j that  utilized  all  n + 1 samples.  Their  estimator  was  developed  for  the 

eciu.il  sample  sizes  case  m,  m,  ...  m ..  Here  the  Korwar  Hollander  est  i mat  or 

11  12  n+ 1 

is  generalized  to  the  unequal  sample  sires  case,  and  exploited  for  simultaneous 


estimation  of  F ,F 


We  establish  a necessary  and  sufficient  condition 


for  the  proposed  estimators  to  be  better  (viskwisc)  than  using,  for  each  distri 
hut  ion,  the  corresponding  empirical  distribution  function  as  an  estimator.  One 


sufficient  condition,  easily  obta  incd  t ruin  our  necessary  and  sufficient  condition, 
is  n * m i n (.  n>  j . •'>  m^  ( j ) > mux  (in  , in , , . . . ,m,,  f j)  • 

2.7.1).  Potential  use  to  Air  Force. 

In  Section  l of  tiiis  paper  we  show  how  the  proposed  estimators  can  be  used 
foi  the  problem  of  simultaneously  estimating  n I distribut  ion  functions.  1 he 
estimators  are  computed  for  some  data  of  Prosehan  (l!)o3)  on  the  times  of  successive 
failures  of  the  air  conditioning  systems  of  Hoeing  720  jet  airplanes.  Of  course, 
the  method  is  applicable  to  survival  analysis  of  other  types  of  systems. 

2. 7. e.  Status  of  paper . 

The  paper  has  appeared  in  The  Theory  and  Applications  of  Reliability,  Vol . I , 

I.  N.  Shiini  and  C.  P.  l’sokos  editors.  Academic  Press,  New  York,  1977. 

2.8.  Testing  whether  more  failures  occur  later. 

2 . 8 .  a . Abst  ract . 

This  is  an  expository  paper  which  illustrates  the  role  of  superadditivity, 
in  reliability  theory,  as  a means  of  describing  wearout  . Statistical  inference- 
procedures  are  described  for  testing  that  minus  the  logarithm  of  a life  distri- 
bution is  superadditive  (i.e.,  the  underlying  life  distribution  is  now  better 
than  used  (NBU)),  and  for  testing  t hat  the  mean  value  function  of  a non- 
homogenous  Poisson  process  is  superadditive.  These  tests  are  appropriate  when 
the  underlying  processes  suggest  that  new  items  are  hotter  than  used  ones  but 
where  we  do  not  insist  that  the  failure  rate  is  increasing. 

i.S.h.  Potential  use  to  Air  force . 

Replacement  policies  are  essential  to  the  maintenance  of  Air  Force  equipment. 

I he  NBU  distributions  play  a fundamental  role  in  the  study  of  replacement  policies. 
Two  replacement  policies  commonly  employed  and  thoroughly  studied  are  age  and 
Mock  replacement.  Under  an  a,  - f , a unit  is  replaced  upon 
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failure  or  upon  reaching  a specified  ago  I,  whichever  route-,  first.  Under  a 
bloak  tvpUia.  -■h,n;  [ a replacement  i •.  made  whenever  a failure  occur.,  and 

add  i t ional  l>  at  specified  rimes  '1,21,51,...  . As  a typical  result,  Marshall 
and  l'rosehan  show  that  a necessary  and  sufficient  condit  ion  for  failure -free 
intervals  to  be  stochastically  larger  (smaller)  under  age  replacement  than  under 
a pol ic\  of  replacement  at  failure  only  is  that  the  underlying  distribution  be 
NB1I  (NKU) . In  a similar  vein,  a necessary  and  sufficient,  condit  ion  that  the 
number  of  failures  in  a specified  interval  [0,  tj  be  stochast i call)  smaller 
(larger)  under  age  replacement  than  under  a policy  »f  replacement  at  failure  only 
is  that  the  underlying  distribution  be  Mill  (NWII)  . Similar  comparisons  bold  using 
block  replacement.  Thus  in  reaching  a decision  as  to  whether  to  use  on  age 
(block)  replacement  policy  or  not,  it  is  important  to  know  whether  the  underlying 
distribution  is  Nlitl  or  not,  equivalently,  whether  or  not  minus  the  logarithm 
of  the  life  distribution  is  superadditive. 

2.S.C.  Status  of  paper. 

This  paper  has  appeared  in  the  Proceedings  ot  the  1978  \pnual  Rcl iahility 
and  Maintainability  Symposium. 

2.9.  Testing  to  determine  the  underlying  distribut ion  using  randomly 
censored  data. 

2. 9.  a.  Abstract . 

For  right  censored  data,  we  develop  a goodness-o f- f i I procedure  for  testing 
whether  the  underlying  distribution  is  a specified  function  Our  tost  statistic 
C is  the  one- sample  limit  of  lvfron's  (I9(>7)  two-sample  statistic  W.  The  tost 
based  on  C is  compared  with  recently  proposed  competitors  due  to  koiol  and 
Green  (1979)  and  Hyde  (1977).  The  comparisons  arc  on  the  basis  of  (i)  applica 
bility,  (ii)  the  extent  to  which  the  censoring  distribution  can  affect  the 


’ 


1 1 
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inference,  and  (iii)  power.  It  is  shown  that  in  certain  situ.it  ions  the  i'  list 
compares  favorably  with  the  tests  of  Koriul -Green  and  Hyde. 

2.9. b.  Pot  on  t i a l use  to  Air  Force. 

The  need  to  generalise  the  goodness -of- fit  problem  to  censored  data  arises 
because  in  some  situations  t lie  observations  are  times  to  the  occurrence  of  an 
end  point  event  and  the  data  are  analysed  before  all  items  have  experienced 
the  event.  For  example,  the  Air  Force  may  be  testing  a given  component  (system, 
etc.)  to  see  if  its  life  distribution  can  be  approximated  by  an  exponential 
distribution  with  mean  1000  days,  and  it  may  be  necessary  to  make  a decision 
before  all  of  the  components  have  fa  lied . 

2.1).  c.  Status  of  paper . 

This  paper  has  boon  submitted  to  Biometrics. 

2. 10.  On  the  normal  convergence  of  a family  of  simple  epidemics. 

2. 10.  a.  Abstract. 

We  consider  a family  of  simple  epidemics  not  necessarily  with  exponential 
inter  infection  times.  For  this  family  of  simple  epidemics  we  establish  normal 
approximations  to  all  finite  joint  state  probabilities.  To  illustrate  the 
applicability  of  our  result  we  present  a class  of  simple  epidemics  used  frequently 
for  modeling  purposes  for  which  the  normal  approximations  hold. 

2.10. b.  Potential  use  to  Air  Force. 

L,et  a filter  in  an  air  conditioning  system  have  N’  orifices  that  can  lie 
clogged  by  particles.  Applying  the  results  of  the  paper  we  can  approximate 
the  number  of  clogged  orifices  at  any  point  in  time. 

2.10. C.  Status  of  paper . 


been  submitted  to  the  dournal  of  Mat hemn t iea I Biosciences. 
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2.11.  The  discrete  asymptotic  behaviour  of  a simple*  batch  epidemic  proce;  >. 

2. it. a.  Abstract . 

A simple  epidemic  process  in  which  the  number  of  individuals  who  can  become 
infected  at  any  point  in  time  is  itself  a random  variable  is  described.  The 
discrete  asymptotic  behavior  of  such  a process  is  disco  sed.  in  particular, 
the  associated  marginal  distribution  and  finite-dimensional  joint  distribution 
of  the  limiting  process  are  considered. 

2.11. b.  Potential  use  to  Air  force. 

Consider  a replacement  of  a tire  (engine,  radio,  etc.)  of  a vehicle. 

Although  the  tire  is  replaced  the  causes  of  failure  will  intensify  over  time. 
Using  the  results  of  the  paper  we  can  approximate  the  number  of  tire  replacements 
up  to  time  t,  even  when  we  allow  more  than  one  source  to  cause  the  failure  of 
the  tire. 

2 . 1 1 .  c . Status  of  paper . 

This  paper  has  been  submitted  to  the  Journal  of  Applied  Probability. 

2.12.  On  the  negative  binomial  convergence  in  a class  of  m-dimensional 
simple  epidemics . 

2. 1 2.  a.  Abstract . 

We  consider  a population  which  is  exposed  to  m infections,  and  consists 
initially  of  N susccpt ibles.  At  each  point  in  time  at  most  one  susceptible 
becomes  infective,  and  only  from  one  cause.  This  m-dimensional  simple  epidemic 

is  a stochastic  process,  (X  (t) X.  (t)),  with  components  counting  the 

number  of  infeetives  from  the  respective  causes  at  time  t.  We  show  that  if 


the  transition  rates  of  cause  1 through  in  at  time  t are  given  by 
m 


"iVi"'"  - S 1 * 1 “l 


0 and  if 


Lim  X.,  . (0)  a b.  t {1,2,...},  then  (X,,  (t),...,X  (tl)  converges  as  N > 

N«  N,i  i 1^,1  N.m 


to 


a random  vector  with  independent  negative  binomial  components. 


2.  I 2 . 1> . 


- K> 

Potential  use  to  Air  Force. 

Let  a filter  in  an  air  conditioning  system  have  N orifice:;  that  can  be 
clogged  by  m different  types  of  particles.  Applying  the  results  of  the  paper  we 
can  approximate  the  number  of  orifices  clogged  by  the  different  types  of  particle 
at  any  point  in  time.  This  application  can  of  course  he  extended  to  other  types 
of  system*  with  ventilation  units,  anti  pollution  devices,  safety  latches,  etc. 

2.12. C.  Statu*  of  paper . 

This  paper  has.  been  submitted  to  the  Journal  of  Stochastic  Processes  and 
t hoi r App l iea t ions. 

2.13.  The  symmetric  m-dimensional  simple  epidemic  process. 

2. 13.  a.  Abstract . 

In  the  classical  simple  epidemic  process  [Bailey  (1975)],  only  one  disease 
is  involved.  In  this  work,  the  model  is  generalised  to  the  case  in  which  sus- 
ceptiblcs  are  simultaneously  exposed  to  m,  m > 1,  diseases.  An  individual  once 
infected  by  a particular  disease  remains  infected  and  is  immune  to  further 
infection  from  another  disease.  The  process  is  called  an  in-dimensional  simple 
epidemic  process.  The  limiting  finite-dimensional  joint  state  probabilities 
and  the  marginal  probabilities  associated  with  each  disease  arc  obtained. 

2.15.1).  Potential  use  to  Air  force. 

Let  a Filter  in  an  air  conditioning  system  have  N orifices  that  can  he 
clogged  be  in  different  types  of  particles.  Applying  the  results  of  the  paper 
we  can  approximate  the  number  of  orifices  clogged  by  the  different  types  of 
particles  at  any  point  in  time.  Again,  this  application  can  be  extended  to 
other  types  of  systems  with  ventilation  units,  safety  latches,  etc. 

2.13. C.  Status  of  paper. 

This  paper  has  been  submitted  to  ’Hie  Royal  Statistical  Society  (Series  B) . 
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2 . i I . Characterizations  of  nonpararaotric  classt  • Lif  dl  stri  but  ions. 

2 . 1 4 . a . Al>  st  ray  t . 

In  this  1 1 1 [ » i»  r we  obtain  charact  or  i -at  ions  of  large  cla  . . of  nonparamutric 

life  distributions,  such  as  the  increasing  (decreasing)  failure  rate,  increasing 
(decreasing)  failure  rate  average,  new  better  (worse)  titan  used,  etc.,  classes. 

The  methoJs  used  differ  from  the  usual  functional  vuu.it  ion  methods  used  for  the 
fnr  more  common  characterizations  of  parametric  families  of  life  distributions. 
2.1't.h.  1’otentinl  use  to  Air  Force. 

The  above-mentioned  classes  of  life  distributions  pin}-  central  roles  in 
the  development  of  tests  and  estimators  for  reliability  data  analysis.  Having 
new  characterizations  of  such  classes  will  provide  new  insights  into  developing 
better  estimators  and  tests  which  can  then  in  turn  be  used  in  Air  Force  appl i - 
cat  ions . 

2.14. C.  Status  of  p.-ijier. 

This  paper  lias  been  submitted  to  the  Annals  of  Probability. 

2.  IS.  On  a charact erization  of  multivariate  distributions  with  applications 

in  reliability  and  epidemiology. 

2 . 15.  a . Ab st  ract . 

let  T T be  positive  random  variables  with  finite  means.  Further  let 

I be  the  collection  of  all  subsets  of  {l,...,n),  and  let  f be  a function  from 
the  nth  Intel  idian  space  to  /.  that  equals  to  J,  (.1  t /)  at  (a  ^ a ) iff 

min  u.T.  < min  a.T..  We  prove  that  min  a.T.  and  (a.,..., a ) are  independent 

. , t i ■ i , i t ...  it  1 n 1 

l cJ  I}:.)  l<isn 

random  variables  for  every  n real  numbers  a,,..., a iff  for  every  n positive 

I n 1 1 

real  numbers  b. and  r l,...,n  the  random  variables  min  a.T. /Ft  rain  a.T.) 

1 n ...  ii  ....  ii 

Isisn  l<i<n 

and  T /F.T  are  identically  distributed.  Further  we  provide  an  explicit  formula 

for  tlu  distribution  of  r.(a, a ).  Multivariate  distributions  that  possess 

I n 1 
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3.  Act ivit ies  ind  honors  of  the  principal  investigator 
during  the  contract  period. 

In  1977  M)les  Hollander  receivoJ  the  honor  oi  becoming  an  elected  member 
of  the  International  Statistical  Institute.  In  April,  1978,  Hollander  was 
elected  by  his  colleagues  to  become  Chairman  of  the  Department  of  Statistics 
at  l-'lorida  State  and  on  July  1,  1978  he  assumed  that  position.  Hollander's 
travel  activities  during  the  contract  period  included  a March,  197'/  \i  it  to 
Madison,  Wisconsin  to  lecture  in  the  colloquia  series.  Department  of  Statistics, 
University  of  Wisconsin;  a May,  1977  visit  to  Columbia,  South  Carolina  to  lector 
in  the  University  of  South  Carolina  Mathematics  Department's  colloquia  series 
and  to  give  an  address  to  the  South  Carol  ilia  Chapter  of  the  American  Statistical 
Association;  a December,  197<>  trip  to  Stanford,  California  to  discuss  potential 
problems  in  censored  data  analysis  with  Byron  Win.  Brown,  Jr.  of  Stanford 
University;  a July,  1977  invited  lecture  at  the  Gordon  Research  Conference 
on  Statistics  in  Chemistry  (New  Hampton,  New  Hampshire);  December,  1977 
invited  lecture  at  the  2nd  International  Symposium  on  Optimisation  Methods 
in  Statistics  (Bombay,  India);  a January,  1978  invited  lecture  at  the  Annual 
Reliability  and  Maintainability  Symposium  (Los  Angeles,  California);  and  an 
August,  1978  contributed  paper  presentation  at  the  lls‘  annual  meeting  of 
the  Institute  of  Mathematical  Statistics. 
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tvorkiiij;  with  Mvlcs  Hollander  on  joint  research  during  t It  contract 
period  were  frank  Proschan,  Jayaram  Sethurnman,  Ranesh  Korwar  ami  C.reyory 
Campbell.  In  the  summer  of  l ‘>78,  Naftali  t anchory,’ s research,  was  supported 
by  this  contract  and  he  undertook  joint  research  effort with  frank  Prosch.in, 
fynne  Hillard,  Rum6n  Icon,  and  Herbert  I.acayo.  Future  joint  research  of 
Myles  Hollander  and  his  current  I’ll. I).  students  Robert  llannum  and  Yuan-Van  dhen 
will  appear  under  AFOSR  78-3678. 


